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We theoretically investigate the transport and magnetotransport properties of three-dimensional 
Weyl semimetals. Using the RPA-Boltzmann transport scattering theory for electrons scattering 
off randomly distributed charged impurities, together with an effective medium theory to average 
over the resulting spatially inhomogeneous carrier density, we smoothly connect our results for the 
minimum conductivity near the Weyl point with known results for the conductivity at high carrier 
density. In the presence of a non-quantizing magnetic field, we predict that for both high and low 
carrier densities, Weyl semimetals show a transition from quadratic magnetoresistance (MR) at low 
magnetic fields to linear MR at high magnetic fields, and that the magnitude of the MR >10 for 
realistic parameters. Our results are in quantitative agreement with recent unexpected experimental 
observations on the mixed-chalcogenide compound TlBiSSe. 
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I. INTRODUCTION 

Electronic band structures that have protected gap¬ 
less points - where the conductance and valence bands 
are guaranteed to meet - have been of significant the¬ 
oretical and experimental interest in recent years. The 
two dimensional manifestation of such band structures 
have been extensively studied in graphene [1], where 
the gapless nature is protected by sublattice symme¬ 
try [2], and in 3D topological insulators [3], where the 
crossing point is protected by topology [4, 5]. More re¬ 
cently, attention has focused on the three dimensional 
analogues of these compounds, called Weyl semimet¬ 
als [6-8]. Compounds such as Cd 3 As 2 [9-11] TaAs [12] 
and TlBiSSe [13, 14] have been shown to have Weyl 
points in their band structure (see Ref. [15] for a recent 
review on the various candidate materials for semimet¬ 
als with 3D relativistic electronic dispersions). 

Theoretical efforts toward characterizing the elec¬ 
tronic properties of Weyl semimetals are in the nascent 
stage and include the scattering properties of differ¬ 
ent impurity potentials [16], localization and delocal¬ 
ization [17, 18], thermoelectric properties [19], screen¬ 
ing [20, 21] and temperature dependence [22], the influ¬ 
ence of the chiral anomaly [23], diffusive transport [24] 
and the effects of electron-electron interactions [25]. 
Inspired by unexpected observations in recent trans¬ 
port [11, 26] and scanning probe [27] experiments, we 
study theoretically the transport and magnetotrans¬ 
port properties of 3D Weyl semimetals in the presence 
of randomly distributed Coulomb impurities. The ef¬ 
fect of the charged impurities is twofold: they provide 
a momentum relaxation mechanism and they act as 
dopants for the local carrier density. At low carrier den¬ 
sity, the former mechanism introduces macroscopic in¬ 
homogeneities in the carrier density profile, giving rise 
to positively or negatively charged puddles. 
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In this work we use a random phase approximation 
(RPA) method to evaluate the effective screened poten¬ 
tial of the Coulomb impurity that enters various disor¬ 
dered averaged quantities. We show that the RPA is 
a much better approximation of the commonly used 
Thomas Fermi approximation due to the nature of the 
vacuum screening structure of Dirac materials. In the 
homogeneous regime (far from the Weyl point in mo¬ 
mentum space), weak impurity scattering is considered 
at the Born level to obtain the Drude conductivity. In 
the inhomogeneous regime (close to the Weyl point in 
momentum space), we find that it is important to con¬ 
sider structure in the disorder distribution when per¬ 
forming the disorder average. We then use an Effective 
Medium Theory (EMT) to average over the inhomoge¬ 
neous carrier density distribution. This formalism has 
been remarkably successful in providing a quantitative 
understanding of the transport properties close to the 
Dirac point in graphene [1], the 2D cousin of these Weyl 
semimetals. Our results allow us to make both quan¬ 
titative and qualitative comparisons with experiment 
and yield many qualitative insights into the behavior 
of Weyl semimetals under various experimental condi¬ 
tions. 


This paper is organized as follows: In Section II, we 
discuss the Drude conductivity using both the Thomas- 
Fermi and RPA screening approximations, far away 
from the Weyl point and also discuss why the RPA 
is required (unlike in the case of graphene). In Sec¬ 
tion III, we discuss the effects of impurity correlations 
and induced charge carriers, that play a role near the 
Weyl point. Finally, we look at experimental results 
from Ref. [26] and compare them with our theoretical 
models in Section IV. The experimental transport data 
is found to be consistent with two possible theoretical 
regimes and we note in Section. V that magnetotrans¬ 
port provides a simple and experimentally accessible 
mechanism to distinguish between the two possibilities. 
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II. TRANSPORT AT HIGH CARRIER 
DENSITY 

The Hamiltonian for a Weyl semimetal is given by 


H = ^hlVF CT - dr — ^ ^V{y) (1) 

Mmp 

V{v)=Y^ C/(r-Rj) 
i=i 


where cr is a vector of Pauli matrices, ^ is the chem¬ 
ical potential, vp is the Fermi velocity and ± ac¬ 
counts for the two chiralities. The density of states 
is i'{E) = where g is the the degeneracy 

(here ^ = 4 due to spin and the presence of two cones). 
In Eq. (1), t/(r-Rj) is the total screened potential seen 
by an electron at position r due to charged impurities 
at positions Hj. In this work, we consider Coulomb 
impurities with momentum space screened potential 


c^(q) 


47re^ 
e(q)q2 ’ 


( 2 ) 


where e is the electronic charge and e(q) is the dielectric 
function. Here k and k' are the incoming and outgoing 
momenta of the scattered electron and q = k — k' is the 
transferred momentum. Note that the Coulomb poten¬ 
tial suppresses large momentum scattering connecting 
the two Weyls cones. For this reason, in the following 
we will work with one cone and consider the contribu¬ 
tion of the second one in the degeneracy factor. For 
a given concentration of impurities, nimp, the ensem¬ 
ble averaged transport scattering time within the Born 
approximation is given by 

h f , 2l-cos26» 

— 27rninip / , .g t^(|k k |) S{E]^ -E-k')- 

Tfr J (27r j 2 

(3) 

To make connection with existing results in the litera¬ 
ture [22, 28] we first consider the simpler case of eval¬ 
uating Eq. (3) in the Thomas-Fermi (TF) approxima¬ 
tion, where the scattering potential can be taken as 
U{q) = + q^)) with a transferred momen¬ 

tum, q = |k — k'|. Here k is the dielectric constant of 
the material and qs = yj^.'Ke^v{EF )/k is the inverse of 
the Thomas-Fermi screening length. Introducing the 
effective fine structure constant a — ^ jhvpK (where 
a = 0.07 for Cd 3 As 2 and a = 0.68 for TlBiSSe), we 
have [22, 29] 


cttf 


e^vjpT 
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^{Ef) 


g kp 1 


(4) 


where H{z) = + 1/2) log(l — 1. In both non- 

chiral two dimensional electron gases and graphene, 
the accidental coincidence at zero temperature of the 
polarization function for q < 2k f and the density of 
states implies that the Thomas-Fermi approximation 
gives identical results to the RPA. However, we should 
emphasize that this is no longer true at finite temper¬ 
ature or in other materials like bilayer graphene [1]. 



FIG. 1. Conductivity as a function of carrier density. 

RPA-Drude theory (Ref. [22], blue curve), self-consistent 
Thomas-Fermi (Ref. [28] , green curve) and the EMT ap¬ 
proach (assuming uncorrelated impurities in both regimes) 
discussed in the main text (red curve). The EMT predicts 
a minimum conductivity close to the Weyl point and repro¬ 
duces the Drude one a ~ far away from the Weyl 

point. In order to compare with the results of Ref. [22], 
we follow their parameters for this figure and use a — 1.2, 
g — 2 and nimp = and liy — 0. (Inset) Close-up 

of the Weyl point. 


As we discuss below, for Weyl semimetals the Thomas- 
Fermi and RPA give qualitatively different results, and 
in what follows, we use the more correct RPA approx¬ 
imation with qs{q) = [^D{E)Ii{q/2kF)Y^‘^ • Here, 
n(x) is the ratio of the RPA polarization function and 
the density of states, and is given by the sum of two 
components: a vacuum part and a finite density 

part [20] 


H 


M 


Hr 


1 + T(1 - 3a;^) log 
4a; 

1 -h X 

-ylog 


1 — X 

x‘^ 

2 

-log 

5 

X 




(5) 


We remark that the Thomas Fermi calculation is 
merely the linearized result of the RPA which is valid 
for q kp [30]. Since we consider transport both near 
and away from the Weyl point, the TF approximation is 
insufficient. Moreover, the TF approximation neglects 
the vacuum term of the polarization function, which in 
our case is the dominant term for g ^ As seen in 
Figure 1, the large quantitative difference between the 
two approximations implies that one must perform the 
full numerical RPA calculation for accurate compari¬ 
son with experiments. Finally, we would like to remark 
that unlike the case of one and two dimensions, for Weyl 
fermions in d = 3 the vacuum polarization function is 
divergent and needs an ultraviolet momentum cutoff A, 
where in Eq (5), d = A/2 /cf. In principle, the transport 
coefficients calculated within the RPA approximation 
could depend on the choice of the cutoff [31], although 
in practice we have verified that such dependence is 
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weak for a realistic parameter range. Our results for 
the high density transport in the absence of a magnetic 
field are shown in Fig. 1. We find that far away from 
the Weyl point, a ^ where n is the carrier 

density. In Fig. 1, we show that this result is in agree¬ 
ment with calculations recently reported in Ref. [22] for 
the homogeneous regime. 


III. TRANSPORT AT LOW CARRIER 
DENSITY 


A. Correlations in Impurity Positions 


In the regime near the Weyl point, we consider two 
different effects. First, the conductivity may be mod¬ 
ified from the homogenous case due to the presence 
of correlations in the impurity positions. In order to 
understand the origin of these correlations, one notes 
that the effective “size” of the impurity potential is ap¬ 
proximately ro — A where A is the screening 

length. This should be compared with the average dis- 
tance between the impurities, given by L 
high carrier density, screening is more effective and the 
screening length is smaller, corresponding to tq ^ L 
(or n ^ Tiimp)- In this regime, impurities are well 
separated and can be essentially considered as point- 
like, meaning that disorder is completely random. In 
the inhomogeneous regime however, we find that the 
electron density is comparable to the impurity den¬ 
sity and therefore rp ^ L. In this case, disorder can¬ 
not be considered as completely structureless and cor¬ 
relations between impurity positions need to be in¬ 
cluded [32, 33]. Since the pair correlation function 
^(Ri,Rj) 7 ^ the self energy term propor¬ 

tional to cannot simply be renormalized away. Fi¬ 
nally, assuming that the impurities are homogeneously 
and isotropically distributed, the pair correlation func¬ 
tion reduces to the radial correlation function g{R). 
The effect of correlations can be taken into account by 
the modified correlator of the Gaussian random fields 
(see Appendix A) 

{V{r)V{r'))d=nimp j U‘^{q) S{q) (6) 

S{q) = 1 + J cfR {g{R) - 1} 

where stands for disorder average and we have in¬ 
troduced the structure factor S{q). Following Ref. [33], 
we take the radial distribution function as 



R< To 
R> ro 


(7) 


When computing transport properties, the effect of 
Eq. (6) is to replace U^{q) U\q) S{q) in Eq. (3). 


B. Effective Medium Theory 

The second effect that one must take into account in 
the inhomogeneous regime is that non-uniformly dis¬ 



FIG. 2. Comparison of n*/nimp obtained using RPA 
and Thomas-Fermi approximations. The error intro¬ 
duced by using the Thomas-Fermi screening to calculate the 
effective minimum carrier density, n*, increases with an in¬ 
crease in the effective fine structure constant a. We note 
that for materials such as Cd 3 As 2 , this is a relatively small 
error but this is not the case for TlBiSSe. Here we use 
cutoff, (5 = 10 


tributed impurities result in a spatially varying local 
chemical potential, which induces a position dependent 
carrier density. The induced carriers in turn screen the 
local potential and one eventually obtains a self con¬ 
sistent relationship between the local induced carrier 
density and the disorder averaged impurity potential, 
V{r,n{r)) = hvF{Q7r‘^n{r)/g)^^^. Note that if we are 
far away from the Weyl point, the fluctuations in car¬ 
rier density are negligible compared to the total num¬ 
ber of carriers but close to the Weyl point this is not 
the case. The average induced carrier density corre¬ 
sponds to the average value of the random Gaussian 
field, Vq = {V{r))d- The net effect of charge doping is 
the appearance of macroscopic regions of charge pud¬ 
dles, each having an excess or a deficit of charge with 
respect to the average value. Therefore, one performs 
an averaging over these spatial fluctuations in carrier 
density using an effective medium theory (EMT). The 
EMT is a well established technique developed by Brug- 
german [34] and later by Landauer [35] in order to char¬ 
acterize the effects of macroscopic fluctuations on the 
global conductivity. In order to perform the average, 
one considers each macroscopic region of definite local 
conductivity to be embedded in a homogenous effec¬ 
tive medium, whose conductivity is determined in a self 
consistent way over all the regions. Early EMT mod¬ 
els assumed two types of regions with conductivities 
a A and as occupying area fractions p and 1 — p. This 
model was later generalized in the case of 2D materials, 
to continuous distributions of local, tensor conductivi¬ 
ties [36, 37]. Here we generalize the results of Ref. [37] 
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to the continuous three dimensional case and obtain 



where is a functional measure, a is the local con¬ 
ductivity and is the effective medium conductivity 
to be found self consistently. In obtaining Eq. (8) we 
have assumed an isotropic material i.e. = ayy = 
and that the local conductivity regions are spheri¬ 
cal in shape [37]. This assumption is valid in the case 
of puddles that are small compared to the sample size. 
The probability distribution in Eq. (8), Vq, K-ms]? 
is the same one that has been used to evaluate the 
Drude transport time. Therefore, it is characterized 
by the average and the variance of the disorder distri¬ 
bution (see Appendix B ). Einally, we define the car¬ 
rier density associated with the variance fd-ms to be n* 
which is the effective minimum carrier density in a Weyl 
semimetal. In Eig. 2, we show that only for a ^ 1 does 
the n* calculated within the RPA reduce to that of the 
Thomas-Eermi. In general, as in the case of transport, 
the full RPA polarizability function must be used. 



FIG. 3. Comparison with the experimental data of 
Ref. [26]. The experimentally determined mobility of Tl- 
BiSSe, /Xexp, decreases as a function of the measured carrier 
density riexp, whereas the Drude theory with constant riimp 
(blue curve) would predict the opposite trend. In the main 
text, we propose two possible scenarios compatible with this 
behavior. The red curve assumes that the experiments are 
in the inhomogeneous regime with riexp ^ rirms with corre¬ 
lated charged impurities. Alternatively, the green curve as¬ 
sumes that the experiments are in the homogeneous regime 
n-exp ^ n-rms, with the charged impurities also responsible 
for doping. Here riexp = no ~ 4 nimp, = 0.68 and (5 = 10. 


C. Crossover Between Homogenous and 
Inhomogeneous Regimes 

In Eig. 1, we show the conductivity at zero mag¬ 
netic field as a function of carrier density. We find that 
for large carrier density, our results reproduce those of 
Ref. [22] but at zero carrier density, and in the inhomo¬ 
geneous regime, our results differ from both the con¬ 
ductivity calculated in Ref. [28] and Ref. [22]. This is a 
consequence of using the EMT which provides a smooth 
crossover between the two regimes instead of a hard 
floor as done in Ref. [28]. Note that Ref. [22] ignores 
the effect of inhomogeneity altogether. The transport 
theory for Weyl fermions shown in Eig. 1 represents the 
full crossover from the inhomogeneous regime (close to 
the Weyl) point to the homogeneous regime (far away 
from the Weyl point). 


IV. COMPARISON WITH EXPERIMENTS 


Next, we apply our theory to address recent experi¬ 
mental findings on TlBiSSe (g = 4, a = 0.68). Ref. [26] 
describe their results of a large increase in mobility with 
decreasing carrier density as surprising. In order to 
make connection with the experimental data, one iden¬ 
tifies the experimental mobility /iexp and carrier density 
Uexp from experimentally determined parameters [38] 


/^exp — 


C7xy(B) 


B-^0 i?(7xx(-B) 


nexp = (10) 


The implicit assumption made here is that the impurity 
concentration does not change with the sample. Indeed, 
for fixed riimp, Fig. 1 shows that mobility increases with 
increasing carrier density. This is represented in Eig. 3 
as the blue curve which shows the opposite trend to the 
experiment. We propose two possible scenarios that al¬ 
low US to relax the constant riimp in a physically justi¬ 
fiable way (and as we discuss below, it is not possible 
to determine from this data alone which of these two 
scenarios correspond to the experimental situation). 

In the first case, we consider the experiment to be in 
the homogeneous regime with charged impurities also 
acting as dopants that shift the chemical potential [29] . 
The average induced carrier density is then given by 
Vo = hvF{Q7r‘^no/gy^^ = nimpt7(q = 0,n*), from which 
we obtain no = 4 nimp in the density regime of inter¬ 
est. The green curve in Eig. 3 uses Uexp = no = 4 nimp 
and shows good agreement with the experimental data. 
Therefore, a plausible resolution of this experimental 
“mystery” is that the charged impurities that are re¬ 
sponsible for scattering carriers are also responsible for 
doping the samples. 

A second possibility is that the samples are in the in¬ 
homogeneous regime where no ^ In this regime, 

one equates the fluctuations in the impurity poten¬ 
tial Vj-ms with the band energy to obtain an effective 
minimum carrier density n*. The red curve in Eig. 3 
considers samples with disorder (nimp varying from 
6.1 X to around 1.9 x 10^^m“^). Using the 

EMT equations (8) and the definition of nexp (10)? we 
see good agreement with the experimental data. Note 


Mexp^ 
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that here we must consider impurity position correla¬ 
tions and take the correlation length ro — 1 / ^ 
as discussed earlier. Both the homogeneous and inho¬ 
mogeneous scenarios are generally consistent with the 
scaling law, /iexp ^ This scaling is a consequence 

of cr ^ n^^^lriimp and n ^ riimp- Note that this also 
presents convincing evidence for using Coulomb impuri¬ 
ties as opposed to neutral or point-like impurities since 
the conductivity scales differently with n for the lat¬ 
ter [ 22 ] . We remark that our analysis is completely free 
of fit parameters and the both the homogeneous and 
inhomogeneous theories are in quantitative agreement 
with the data (the inhomogeneous curve agrees with the 
data to within a factor of 4). This also demonstrates 
that transport measurements alone cannot distinguish 
between the two regimes since both show comparable 
mobilities (same order of magnitude). We propose that 
magnetoresistance is the appropriate measurement to 
distinguish between the two regimes. 


V. MAGNETORESISTANCE 

To calculate the magnetotransport of Weyl semimet¬ 
als, we assume that the charge and the axial current are 
weakly coupled, since Coulomb impurities suppress the 
scattering between the two Weyl nodes [39]. Indeed, the 
results of Ref. [26] show no negative magnetoresistance 
or any effect of an in plane magnetic field, both of which 
are signatures of the chiral anomaly [40, 41]. Moreover, 
we note that in the absence of parallel electric and mag¬ 
netic fields, the effects of the chiral anomaly are absent. 
Thus, the origin of transverse magnetoresistance in this 
system is likely due to be disorder-induced which can 
be treated in a semiclassical regime. We then assume 
that one can define a conductivity matrix with 

1 jiB 

= = ( 11 ) 

where ( 75 ( 1 *) = cr(n(r)) is the RPA-Boltzmann conduc¬ 
tivity discussed earlier and the magnetic field is taken 
along the z axis. The above set of equations con¬ 
stitute the input for the EMT model. The solution 
of the EMT equations contains both the magnetore¬ 
sistance caused by having two types of carriers (elec¬ 
trons and holes) [42] and the disorder-induced mag¬ 
netoresistance discussed in the context of silver chalco- 
genides [37, 43] and other two dimensional systems [44], 
where MR = {pxx{B) - pxx{B = 0))/pxx{B = 0). 

A simple physical picture for this disorder induced 
MR was proposed in Ref. [44] and this still holds in 
Weyl semimetals. Essentially, disorder results in a sit¬ 
uation where the charge carriers move with varying 
drift velocities in various regions. The global Hall field 
thus cannot cancel the velocity dependent Lorentz force 
(as would be the case in a homogenous system) and 
the electron trajectories become longer as the magnetic 
field is increased. This is the origin of disorder induced 
magnetoresistance. 



FIG. 4. Magnetoresistance in TlBiSSe. Disordered 3D 
Weyl fermions can have large MR > 10 for realistic ex¬ 
perimental parameters. The figure shows that both in the 
homogeneous and in the inhomogeneous regime, the magne¬ 
toresistance is quadratic at low fields (see inset) and linear 
at high fields in agreement with experimental observations. 
Notice that the MR in the inhomogeneous regime is much 
larger than that of the homogeneous regime, suggesting that 
increasing disorder is an easy way to enhance the MR. 


Eor Weyl semimetals in general, our theory predicts 
that the MR should be quadratic at low magnetic 
fields and linear at high magnetic fields. In Eig. 4 
we show our results for uq = 3.8 x and 

^imp = 9.5 X (homogenous regime), and 

for no ^ 0, nimp = 2.3 x (inhomogeneous 

regime). These values were chosen so that they 
correspond to similar /iexp and Uexp and therefore they 
cannot be distinguished from transport measurements 
alone. The results in Eig. 4 demonstrate convincingly 
that within the semi-classical theory presented here, 
the magnetoresistance in the inhomogeneous regime 
is much larger than that of the homogeneous regime. 
We also note that in the inhomogeneous regime MR 
is comparable to that seen by Ref. [26] for matching 
parameters, suggesting that those samples were in the 
inhomogeneous regime. Within this model, having 
MR > 10 is easily achievable in the inhomogeneous 
regime for moderate values of B, but it is several orders 
of magnitude weaker in the homogeneous regime. This 
suggests a clear way to experimentally distinguish 
between these two regimes. Moreover, an easy way 
to increase MR for technological applications is to 
make the sample dirtier, a counterintuitive, yet easily 
achievable, experimental goal. 
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X ^ ^ 

s=<i+-L_L>=i. + /' \ + i i 

FIG. 5. Disorder averaged self energy. Amputated di¬ 
agrams defining the self energy at the Gaussian level. The 
brackets represent disorder averaging, the “x” accounts for 
the impurity concentration, the dot for the impurity poten¬ 
tial insertion and the solid line is the clean propagator. 

their experimental data. 

Appendix A: Correlated Impurities 

Here we consider the effect of spatial correlations be¬ 
tween impurities in the inhomogeneous regime. This 
problem has been considered before in Refs. [33, 45], 
but its relation to diagrammatic perturbation theory is 
now made explicit. Disordered averaged Green’s func¬ 
tions are defined in terms of their self energy. At the 
Gaussian level, the disorder averaged self energy is de¬ 
fined in terms of the diagrams depicted in Fig. 5. Disor¬ 
der averaging is generally defined in terms of a weighted 
sum over impurities position of a disorder dependent 
quantity 0(Ri) [46] 

V %mp 

n dR*5(Ri)0(r,Ri) (Al) 

i = l 

N- N- 

^ zmp ^ ^ imp 

n dKi dKj gCRi.Kj) (A2) 

i=i j=i 

X 0(r,Ri)0(r',R^) 

Here (...) stands for disorder average, R^ is the position 
of the impurities in a d-dimensional volume Ni^p 

is the total number of impurities in the system and the 
indices i, j, /... label different impurities. The important 
objects in the above definitions are the correlation func¬ 
tions ^(Ri, Rj,..., R;z), describing correlations between 
one impurity, two impurities and so on [32]. As usual, 
the hierarchy of correlation functions cannot be worked 
out explicitly and one has to perform some physically 
motivated ansatz in order to truncate the hierarchy. 
For completely random disorder, all correlation func¬ 
tions factorize as product of single particle correlations 
[46]. These are simply equal to the probability 
of finding an impurity at site i, i.e. 1/L^. The pair 
correlation function ^(Ri,Rj) gives the probability of 
finding an impurity at site Ri given one at site Rj . Let 
us consider the total disorder potential V (r) defined in 
Eq. (1) 

l/(r) =f^U{r- R,). (A3) 

i=l 

In the weak scattering limit (Born limit), only the first 
two moments of the distribution of E(r) are relevant. 


corresponding to the diagrams of Fig. 5. Performing 
the impurity average, one obtains three diagrams: two 
proportional to Ni^p and one to For completely 

random disorder, the term proportional to is a 

constant that can be renormalized away. However, 
as we are now going to show, if there is any residual 
structure in the disorder, the term cannot be 

renormalized away. 

For completeness, we start considering the standard 
single impurity diagram; in momentum space it reads: 

N- 

^ ^ imp 

^ (C/(k - U{0) <5(k - k'), 

i=l 

(A4) 

where riimp = Ni^pjL^ is the density of impurities and 
U (0) is the impurity potential evaluated at zero trans¬ 
ferred momentum q. Note that in principle this term 
is singular and needs to be regularized, for example 
by screening. Here we will assume that the impurity 
potential is screened and therefore consider U{0) as a 
finite quantity. Consider now the second order term 

^ imp ^ 

i,j=i k' 

X U{k' - k") e-*R-Ak'-k")^^ (A5) 

where Go(k', E) is the free electron propagator. Taking 
the disorder average, there are two contributions: for 
i = j one finds 

^ imp 

=ni„p5(k-k"), 

i=l 

(A6) 

corresponding to the rainbow diagram shown in Fig. 5. 
For i jf, we need to make some assumptions on the 
pair correlation function. We assume that it still de¬ 
pends only on the coordinate difference R^ — Rj [32]. 
As a consequence, momentum is conserved on average, 
i.e. k = k". We also assume that the pair correlation 
function does not depend on the angle between R^ and 
Rj. Within these assumptions, Eq. (A5) reads 

^J2(^o{k',E)uWk-k'\)nimp (A7) 

k' 

X f^l + riimp J 

where the first term in the curly bracket corresponds 
to the Born term (single particle scattering) and the 
second takes into account the effect of two particles 
scattering. Note that the second term is singular at k = 
k' [32] ; to take care of this singularity, one subtracts it 
as the Eourier transform of unity and we can define the 
regularised structure factor as 

S'(lk-k'l) = l + Uimp f d‘^R[g{R) - 1] 

(A8) 
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This has been obtained in Refs. [33, 45] and the struc¬ 
ture factor can be measured e.g. in neutron diffraction 
experiments. Next, we consider the case where the term 
proportional to is important. For the charge carri¬ 
ers, the “size” of the charged impurity is roughly given 
by the effective Bohr radius uq of its lowest impurity 
level; this can be orders of magnitude larger [47] than 
the underlying lattice constant [32]. If one compares 
the average spacing L between the impurities and ao, 
one comes at the conclusion that correlation effects are 
important if L < ao- Note that the ratio L/a^ is rem¬ 
iniscent of the parameter used to quantify correla¬ 
tions in an electron gas. We conclude by connecting the 
above analysis to the relaxation time r. By definition 
l/r(/c) = ImS(/c), where 

Inil](fc) = [/2(|k-k'|)5(|k-k'|)ImGo(k',i;) 

k' 

(A9) 

It follows that the variance of the random Gaussian 
field V (x) can be written as 

(V(r)V(r')}d = U\q)S{q) 

(AlO) 

that is our Eq. (6) of the main text. 

Appendix B: On the Gaussian approximation 

Here we consider the form of the disorder probabil¬ 
ity distribution used in Eq. (8) of the main text. As 
explained in the main text, within the Drude trans¬ 
port theory, this function is the same as the one used 
to evaluate the transport time. Here we provide an 
explicit proof of the validity of the Gaussian approx¬ 
imation. This discussion is mostly based on Ref. [48] 
and is based on the functional approach to disordered 
systems. This method is completely equivalent to the 
“sum over impurities” approach used in diagrammatic 
perturbation theory. In the functional approach, one is 


interested in evaluating the disordered averaged gener¬ 
ating functional 

{logZ[V])d = I ^VP[V] logZ[V], (Bl) 

where Q)V is a functional measure. Here we are not in¬ 
terested in the actual calculation of this quantity, but 
only in finding P[E]. Generally, log Z\y\ can be substi¬ 
tuted with an effective functional of V as in the case of 
the Effective Medium Theory (EMT) employed in the 
main text. Let ns consider the total disorder potential 
of Eq. (A3), where the potential U is due to the Nimp 
impurity potentials and E(r) is the resulting effective 
potential. In order to find P[E], one needs to interpret 
the above definition as a constraint. In the functional 
formalism this is accomplished by means of a functional 
Dirac delta function averaged over impurity positions 

^imp 

P[V] = n(<^[l^(r) - ^ f/(r - Ri)]) (B2) 

r i=l 

N • 
imp 

zfPriir)[V{r)- E U{r-Ri)] 

e ), 

where we have used the functional representation of the 
Dirac delta function and ^(r) is a Lagrange multiplier 
field. The first step in the evaluation of P[V] consists 
in assessing the correlated nature of the impurities. Ac¬ 
cording to Ref. [46], and as explained in the main text, 
if the distance between the impurities is much larger 
than the screening length, then the positions of the im¬ 
purities are uncorrelated. On the other hand, if the 
screening length of the Coulomb potential is compa¬ 
rable to the average distance between the impurities, 
then there may be correlations in the positions of the 
impurities themselves (see Appendix A). Here we con¬ 
sider for simplicity the case of completely uncorrelated 
disorder and comment on the effect of correlations at 
the end. In the thermodynamic limit one obtains [48] 



P[V] = 


J- 


fPr^{r)V(r) 


exp < —n 


f^imp 


j d^R(^ 


1 — e 


-^fPr^ir)Uir-R)\ I ^ 


/ 


= / 


J d^r ar)Vir) <^>(0 


(B3) 


where in the second step we have identified the cumu- 
lant function of the stochastic process 4>(^) and the 
related characteristic function xiO — • We now 

follow Ref. [49] (where more details can be found) and 
assess the Gaussian nature of the characteristic func¬ 
tion using the central limit theorem. One needs to ex¬ 
pand the cumulant function and show that the magni¬ 
tude of the second cumulant is the dominant one. This 
is basically the well known Born criterion, c.f. Ref. [46]. 


Eor screened Coulomb disorder, and for the considered 
values of the electromagnetic coupling a, the Born cri¬ 
terion is satisfied for carrier density n > n*, see Sec¬ 
tion HIB. Einally, integrating out the Lagrange mul¬ 
tiplier field, one finds the Gaussian distribution of the 
random fields used in the main text. When n n*, 
the screening length of the Coulomb potential is com¬ 
parable to the average distance between the impurities 
and one needs to use the modified propagator of the 
random fields, Eq. (AlO). 
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